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We study quantum effects of light propagation through an extended absorbing system of two- 
level atoms placed within a frequency gap medium (FGM) . Apart from ordinary solitons and single 
particle impurity band states, the many-particle spectrum of the system contains massive pairs of 
confined gap excitations and their bound complexes - gap solitons. In addition, "composite" solitons 
are predicted as bound states of ordinary and gap solitons. Quantum gap and composite solitons 
propagate without dissipation, and should be associated with self-induced transparency pulses in a 
FGM. 

PACS numbers: 42.50.Md, 42.50.Ct, 42.65. Tg 

The self-induced transparency (SIT) pulses, predicted and observed in the pioneering work of McCall and Halm 
jjj , may be regarded as solitons of the Maxwell-Bloch model , describing classical radiation propagating in a single 
direction and coupled to an extended system of two-level atoms. The model is completely integrable [|| and the time 
evolution of an arbitrary radiation incident on an atomic system is described by the inverse scattering method 

§• 

In the case of a high intensity pulse in ordinary vacuum, quantum corrections are negligibly small. Therefore 
the quantum version of the classical model - quantum Maxwell-Bloch (QMB) model - has been studied [Q only in 
the context of the superfluorescence phenomenon where quantum effects play a crucial role |ql. But the situation is 
drastically changed for frequency gap media (FGM), such as a frequency dispersive medium |9|, a photonic bandgap 
material [jTc| , [TT[ ] , and a one-dimensional Bragg reflector [fl2[ , where classical, linear wave propagation inside a frequency 
gap is excluded @0. 

In this Letter, we demonstrate the existence of nonclassical light propagation through an extended homogeneous 
fL5| system of two-level atoms placed within a FGM. These light pulses are highly correlated quantum many-body 
states and are distinct from single photon hopping conductivity through the photonic impurity band created by 
the atoms. Because of a nonlocal polariton-atom coupling, an extension of the Bethe ansatz method |l7]] from the case 
of a single atom embedded in FGM Jl3|,[l4j to the case of an extended many-atom system requires a thorough analysis. 
The QMB generalized to the case of FGM exhibits hidden integrability []l3| , provided that the characteristic times of 
interatomic resonance dipole-dipole interaction (RDDI) and other collisional dephasing effects are much longer than 
the light pulse duration. Integrability of the model allows us to describe the time evolution of an arbitrary light 
pulse incident on the system in terms of the allowed soliton modes. Making use of the Bethe ansatz technique, we 
derive the Bethe ansatz equations (BAE), which completely determine the spectrum of the radiation plus medium 
plus atoms system. Here we consider the case of the atomic transition frequency u}\2 lying deep inside a frequency 
gap of a frequency dispersive medium, for which the McCall-Hahn theory is inapplicable. 

Unlike an attractive effective photon-photon coupling in empty space (or nondispersive media) caused by scat- 
tering of photons on an atomic system, an effective polariton-polariton coupling in FGM is shown to be attractive 
only for polaritons of the lower polariton branch. Therefore bound many-polariton complexes (ordinary solitons) can 
be constructed only from polaritons of the lower branch. In the limit of a macroscopically large number of polaritons, 
these quantum complexes arc nothing but SIT pulses (27r-pulses) of the classical theory slightly modified due to a 
nonlinear polariton dispersion. 

Due to the existence of a frequency gap, the multiparticle spectrum of the system, apart from polaritons and 
ordinary solitons, also contains massive pairs of confined gap excitations, which do not exist out of pairs, and bound 
complexes of these pairs - quantum gap solitons. The energy-momentum dispersion relations for gap solitons are 
derived and the widths of soliton bands, the soliton masses, the spatial sizes and the velocities of propagation inside 
the atomic system are evaluated as functions of the number of pairs and the atomic density. We predict also the 
existence of "composite" solitons as bound states of "deformed" ordinary and gap solitons. 

In contrast to quantum gap solitons generated by a single atom, which propagate along a radial coordinate centered 
at the atom E3], quantum SIT pulses in a doped FGM propagate in a direction defined by a single wave vector. 
Furthermore, the gap SIT pulse, consisting of an even number of gap excitations, is distinct from (odd photon 
number) gap soliton hopping conduction inside the RDDI mediated impurity band. 

In the dipole, rotating wave approximation || the Hamiltonian of the generalized QMB model can be written as 
H = Hg + V, where 
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represents the Hamiltonians of M identical two-level atoms and free polaritons, while the operator 

v = -V^E / ^h + pMe ,tMl " + h-c] (ib) 

0=1 ^ c 

describes their coupling. The polariton operators p(w) obey the commutator [p(tu),p\tu')\ = 2t:5(lo — w'), while the 
spin operators a a — (a%, a v a , a*), a ± = <r x ± i<r a describe atoms having the coordinates {x a , a = 1, . . . , M} on the 
polariton propagation axis (the x axis). The states between frequencies labeled as f2j_ and are forbidden for linear 
propagating polariton modes. Therefore, the integration contour C consists of two allowed intervals, C = C_ + C+, 
where C_ = (0, Ox) and C+ = (fin, 00). The coupling constant 7 = 2ttloi2<1 2 / 'So, where d is the atomic dipole 
moment and So is the cross-section of a light beam. The information about the medium spectrum is contained 
in the dispersion relation k(uj) = um(u>). The atomic form factor z(uj) = luti 3 (uj)/lo 12, where n{ui) — y/e(w) and 
e{to) = (oj 2 — £l 2 )/(ui 2 — f2j_) is the dielectric permeability of a frequency dispersive medium. 
The eigenvalues of the model (1) are found from the following Bethe ansatz equations 

where kj = k{ujj), E — u>j is the eigenenergy and the "rapidity" hj = h(cjj) is given by 

Kv) - (2b) 

BAE have a clear physical meaning: the first phase factor on 1. h. s. is acquired by a polariton wave function during 
free propagation between points =fL/2, while the second one accounts for phase factors resulting from subsequent 
scattering on M atoms. Propagating between the points the polariton is also scattered by the other N — 1 

polaritons, and its wave function acquires the phase factor given on r. h. s. in eq. (2a). Information concerning 
the nonlinear polariton dispersion is contained in the rapidity h(u>). In empty space, it is reasonable to neglect 
the resonance dipole-dipole interaction between atoms in self-induced transparency, since the dominant photon-atom 
interaction is pointlike and the excited atom decays by stimulated emission into optical pulse modes. In this case, light 
scattering from each of the M atoms is considered independently. In the FGM, there are no classical modes available 
for stimulated emission and the polariton-atom coupling is highly nonlocal: An excited atom (photon-atom bound 
state) exhibits nonlocal interaction with other atoms which are within the classical tunneling distance. This leads 
to coherent hopping conduction of a photon through the resulting impurity band [ ft6[ . If the characteristic time of 
the interatomic hopping is much longer than the pulse duration, RDDI mediated transfer of energy between impurity 
atoms in an arbitrary direction can be neglected. Energy transfer occurs through the soliton band |14j rather than 
the impurity band |lq| . Eq. (2a) is obtained by including RDDI contributions only from virtual polaritons traveling 
in the same direction as the incident pulse. In this generalization of the single-atom soliton band to an M-atom 
soliton band, scattering from each atom is treated independently. This is equivalent to the independent atom (gas) 
approximation used by McCall and Hahn in ordinary vacuum. Accordingly, in ordinary vacuum [n(u>) = 1], eqs. (2) 
reduce to the BAE of the QMB model which, for large TV, describes the classical McCall-Hahn solution. 

It is instructive to derive the main results of standard SIT theory from eqs. (2) with the rapidity h(uj) ~ (lu— ujyi) j 
corresponding to the case of empty space. As L — > 00, eqs. (2) admit solutions in which complex rapidities hj are 
grouped into "strings" containing n particles, 

h j =H + i((3/2)(n + l-2j), i = l,...,n, (3) 

where if is a common real part ( "carrying" rapidity) . Due to the linear relationships between the rapidity, frequency 
and momentum, particle frequencies and momenta are also grouped into string structures, kj = K + i("//2)(n+l — 2j), 
ujj = Q + 2(7/2) (n+ 1 — 2j), where K and O are common real parts of momenta and frequencies, respectively. To avoid 
confusion, we use the term "string" for solutions of BAE in the ft,-space and the term "soliton" to refer to string's 
images in the lu- and fc-spaces. Consider for simplicity the case when all N particles are grouped into a string, i. e. 
N = n. Inserting kj and LOj in eq. (2a) and evaluating the product over j = 1, . . . , N, we obtain the simple equation 
exp (iQnL) = 1, where 
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Q(Q) = K-— arctan — --^ (4a) 

n 2(11 ~ UJ12) 

and the number of atoms is represented as M = pL. Here p is the linear density of the number of atoms. Clearly 
Q(il) can be interpreted as the energy- momentum dispersion relation of a soliton of size n, where the second term 
describes a contribution of photon-atom scattering. The group velocity of soliton propagation V = dQ/dQ is then 
given by 

J_ = 1 IP Uh] 

V c (n-^ 12 ) 2 + (n 7 /2) 2 ' 1 ' 

where c is the speed of light in empty space. Eq. (4b) is identical to the corresponding expression in the classical SIT 
theory. The spatial size of a soliton, l s ~ (771) -1 , is inversely proportional to the number of photons 0. Therefore, 
only macroscopically "long" strings, b> 1, propagate in an absorbing atomic system without dissipation. 

In FGM, eq. (3) is a solution of BAE if and only if the imaginary parts of rapidities hj and corresponding momenta 
kj have the same sign, 

sgn^m/ij) = sgn(Imfcj), j = l ) ...,n. (5) 

It is easy to understand that the necessary condition (NC) (5) determines the frequency intervals, in which an effective 
particle-particle coupling is attractive, and hence admits bound many-particle complexes. 

We start with the case when the real part of ojj lies outside the gap. Let oj = £ + irj and £ € C. Making use of the 
approach developed in jl4j , it is easy to show that the effective coupling is attractive only between polaritons of the 
lower branch, £ € C_. Polaritons of the upper branch are described by one-particle strings with real positive rapidities 
and do not form any bound complexes. Bound many-polariton complexes (ordinary solitons) are quite similar to 
solitons of the QMB model, despite their inordinate behavior on different polariton branches. The dispersion relation 
of an ordinary soliton of size n is given by 

?(£) = fc(f) - — arctan J^rr, (6a) 

where fc(£) = The group velocities inside, V = d^/dq, and outside, v = d^/dk, the atomic system are then 

related by 

1 = 1 + el d m (6b) 

V v h 2 ^) + {l3n/2f d£ 1 ' 

Since ordinary solitons in the FGM are off-resonance to the atomic transition, the effect of the atomic system on their 
propagation is always weak, unlike the case of SIT in empty space. 

Next we study the multiparticle excitations of the system with eigenenergies lying inside the frequency gap. We 
look for an image of a Bethe string when the real parts of particle frequencies u)j lies inside the gap, £ e G = (f2j_,0|i). 
To find the analytical continuations of the functions k(u>) and h{cu), an appropriate branch of the function n(u>) is 
fixed by the condition n(£±iO) = ±w(£), where u(^) = y/\e(£)\. Soliton parameters £j and r)j are determined by the 
equations 

RehfaVj) = H, hnhfe,Vj) = N + 1/2- i). (7) 

Approximate solutions of eqs. (7) are easily found by the method of Ref. One can show that in the model under 
consideration gap excitations exist only for £ € (^12, ^||)- A string with an even number of particles, n = 21, describes 
a bound complex of / pairs of confined gap excitations - quantum gap soliton - with the eigenenergy per pair of 

e^e^ -A l + q 2 /2mi, (8) 

where q is the soliton momentum per pair, while A; = ^r(4^ 2 — 1) is the band width, and m i — % [j Ylj \ K '(£j ) 

is the mass of a gap soliton containing I pairs. Here k' = dn/d£, k(£) = £^(£)> an d Q — Wi2 + (j3/a)(l + l/2 — j). The 
parameters a and b are found as coefficients in the Taylor series of the function 0(£) = [£^ 3 (£)] _1 at the point £ = cj 12 : 
</>(£) ~ a + b(£ — LU12). The positions of the centers of soliton bands are given by = l^ 1 J2j £j — ^12 + (P/ a )l- 
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The spatial size of a pair, 6 ~ k _1 (£), is nothing but the penetration length of the radiation with the frequency 
u> = £ G G into the medium, and hence it lies on the scale of a few wavelengths. Since is monotonically decreasing 
function, the gap soliton size, in a sharp contrast to the case of ordinary solitons, grows with the growth of the number 
of pairs I. 

For a gap soliton, the dispersion relation inside the atomic system takes the form 



The group velocities of gap soliton propagation inside, Vi = dei/dQ, and outside, vi = de/dq, the atomic system are 
then related by 



Again in a sharp contrast to the case of ordinary solitons, the velocity of gap soliton propagation inside the atomic 
system is greater than vi, i. e., the particle-atom scattering speeds up a gap soliton. But, it should be emphasized 
that the results obtained are valid only for quite small soliton momenta q when the gap soliton dispersion is described 
in the effective mass approximation. At arbitrary q or for quite big solitons (I 3> 1), we cannot use only the first 
terms of the Taylor expansion for the function </>(£) and have to solve the exact equations (7). 

So far, we have looked for a soliton image of a string assuming that all of the string rapidities are mapped to the 
soliton frequencies whose real parts lie either outside (ordinary soliton) or inside the gap (gap soliton). It is easy 
to see that this assumption is not necessary, and one can construct bound many-particle complexes - "composite" 
solitons - containing excitations from different frequency intervals C_ and G. To clarify this point, let us consider first 
the simplest example of a three-particle string with a negative carrying rapidity H. Then, the pair of the complex 
conjugated rapidities hi = H + i(3 and h 3 = h* = H — i[3 is mapped to the pair of frequencies u) = £ + it] and 
lu* = £ — ir], where £ = W12 + (3/a G G and r\ = \H\/a, and momenta k = q + i/c(£) and k* = q — where 
q = t)\k'(£)\. Under the condition H < 0, the real rapidity h 2 = H can be mapped to the real frequency G C_ and 
the corresponding real momentum fc_ = where is determined by /i(£_) = H. All three particles compose a 

bound state of a pair of gap excitations and a polariton of the lower branch. 

The generalization of the composite soliton construction to the many-particle case is straightforward. Consider an 
Af-particle string with an negative H. A pair of the complex conjugated rapidities or any number of such pairs could 
be mapped to corresponding pairs of frequencies whose real parts lie inside the gap, while the rest rapidities of the 
string are mapped to frequencies whose real parts lie below the gap. In other words, one part of the string rapidities 
is used to construct a "gap soliton", while the rest rapidities of the same string are used to construct an "ordinary 
soliton" . All the particles together compose a bound many-particle complex, so that an composite soliton can be 
treated as a bound state of an "ordinary" and a "gap soliton" , but these "solitons" are different from ordinary and 
gap solitons discussed above and do not exist separately from each other. 

Both gap and composite solitons realize the self-induced transparency effect in FGM. These solitons describe an 
entirely new mechanism for energy transfer in a frequency gap medium. The results obtained illustrate the rich variety 
of nonclassical optical propagation effects within a classically forbidden frequency gap. Furthermore, they suggest 
that a doped FGM with suitable optical pumping may act as a source of novel quantum correlated states of light. 
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